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Abstract. For the treatment of interacting electrons in crystal lattices approxima- 
tions based on the picture of effective sites, coupled in a self-consistent fashion, have 
proven very useful. Particularly in the presence of strong local correlations, a local ap- 
proach to the problem, combining a powerful method for the short ranged interactions 
with the lattice propagation part of the dynamics, determines the quality of results to a 
large extent. For a considerable time the non crossing approximation (NCA) in direct 
perturbation theory, an approach originally developed by Keiter for the Anderson im- 
purity model, built a standard for the description of the local dynamics of interacting 
electrons. In the last couple of years exact methods like the numerical renormalization 
group (NRG) as pioneered by Wilson, have surpassed this approximation as regard- 
ing the description of the low energy regime. We present an improved approximation 
level of direct perturbation theory for finite Coulomb repulsion J7, the crossing ap- 
proximation one (CA1) and discuss its connections with other generalizations of NCA. 
CA1 incorporates all processes up to fourth order in the hybridization strength V in a 
self-consistent skeleton expansion, retaining the full energy dependence of the vertex 
functions. We reconstruct the local approach to the lattice problem from the point of 
view of cumulant perturbation theory in a very general way and discuss the proper use 
of impurity solvers for this purpose. Their reliability can be tested in applications to 
e.g. the Hubbard model and the Anderson-lattice model. We point out shortcomings 
of existing impurity solvers and improvements gained with CA1 in this context. 
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1. Introduction 

In a key paper [1], H. Keiter and J. C. Kimball in 1971 described a new perturbational 
method for treating the problem of an impurity with strong local Coulomb matrix 
elements, embedded in a metallic host. Their guideline was to preserve the local 
correlations from the outset, contrary to Hartree-Fock theory or to other decoupling 
schemes, and to keep the interpretation of individual contributions as physical processes. 
The particular difficulties to be surmounted arose from the fact, that choosing 
hybridization or intersite transfer of single particles as the perturbation leaves an 
interacting local shell as the unperturbed part of the Hamiltonian. In such a case the well 
known machinery of Feynman diagrammatics cannot be used, including Wicks theorem 
and linked cluster expansions. The solution found used time-ordered pieces of Feynman 
processes, visualized as Goldstone diagrams, and organized them in the form of Brillouin- 
Wigner perturbation theory with real energy variables. In this early formulation of 
the theory the need to regularize vanishing energy denominators prevented extensive 
studies to infinite perturbational orders, which are necessary in the presence of infrared 
divergencies, encountered e.g. in the Kondo problem. These can be thought of as to 
arise from degeneracies in a classical part of the Hamiltonian; they are then lifted 
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by quantum fluctuations. Early applications of the technique can be found in [2], 
where leading logarithmically divergent terms are summed to all orders to generate 
finite results, and in [3, 4], where also generalizations of the formalism to more general 
local shell structures and to the genuine lattice problem were discussed. 

The technique of Keiter and Kimball originally was designed in connection with the 
Kondo problem. A revival occurred around 1980, when metallic compounds exhibiting 
the intermediate valence phenomenon stayed in the focus of experimentalists [5], and 
somewhat later, when the existence of very heavy quasiparticles (Heavy Fermions) 
was revealed in such compounds containing ions with active 4f- or 5f-shells [6]. In 
particular the discovery of Heavy Fermion superconductors [7] spurred the investigation 
of lattice models with strong local correlations. A breakthrough in theory occurred in 
1983, when it was learnt how to handle direct perturbation theory completely in the 
complex plane thus circumventing the regularization problem. The first independent 
extensive studies [8, 9] of impurity problems were based on the leading skeleton 
diagram contributions to the dynamics of ionic shell states, requiring the solution 
of self consistently coupled singular integral equations. The contributing diagrams 
are completely characterized as not containing crossing band electron lines. This 
approximation (NCA) was shown to adequately describe dynamical properties like 
excitation spectra in addition to thermodynamic quantities, the latter of which being 
known from Wilson's implementation of the numerical renormalization group [10]. 
NCA allows to study [8] the temperature dependent formation of the Abrikosov- 
Suhl resonance (ASR) [11] with an assessment of the systematic shortcomings via a 
comparison to the resonant level model, the limit of vanishing spin degeneracy. NCA 
can also be characterized as a conserving approximation in the sense of Baym and 
Kadanoff [12] and can consistently be applied to calculate various properties of the 
strongly correlated impurity problem [9]. Also in 1983, a reformulation of perturbation 
theory was presented [13], which allowed for the use of Feynman diagrammatics via 
the introduction of auxiliary particles (Slave Bosons). Correlated local states are 
reintroduced in this approach by a constraint on the larger Hilbert space and a 
corresponding projection onto the physical sector after resummations. The Slave Boson 
method stands in one to one correspondence to the formulation via direct perturbation 
theory and thus contributes to the same line of development. 

Merits and shortcomings of the NCA for the Anderson impurity with infinite 
local Coulomb repulsion meanwhile are well known, for example as result of the early 
numerical studies [8, 14] or of the exact analysis of the case with a flat electron-hole 
symmetric conduction band at zero temperature [15]. NCA e.g. captures the exponential 
part of the Kondo temperature T K) the dynamically generated energy scale below which 
a local Fermi liquid is formed due to spin compensation, but not the prefactor. It does 
not furnish the correct values of the threshold exponents, connected with the time 
development of ionic states as known from the X-ray absorption problem [16, 17]. Its 
accuracy increases with increasing (orbital) degeneracy of the ionic level, and NCA may 
even become a fully acceptable approximation for a multi-channel situation [18]. 
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The need for improved treatments of the impurity problem with correlated 
electrons, however, turned out to be even more important, when it became a building 
block in theories of lattice systems [19, 20, 21], using the concept of effective sites. It 
turned out that e.g. the coherence forming in the low temperature regime of lattices 
can only be correctly retrieved by a proper incorporation of the effective local Fermi 
liquid. The necessary increase in the quality of results for lattice models was largely 
driven by the use of numerically exact methods as impurity solvers, i.e. quantum Monte 
Carlo (QMC) and numerical renormalization group (NRG) [22]. NRG in particular, 
in its extension of the original static version to dynamical quantities [23], has been 
developed to a useful and convenient tool for a description of the low energy region [24] . 
When combined with the (cluster-) dynamical mean field theory (DMFT) it opens the 
perspective for a proper description of ground states and correlations in the lattice. 

Still there remain problems: (1) Non local matrix elements of the Coulomb 
interaction can be strong, too, and may enforce a nonlocal approach with links 
between sites the outset, as e.g. in pyrochlore lattices, where for certain fillings classical 
ground states characterized by the tetrahedron rule build the arena for quantum 
fluctuations [25]. (2) The use of clusters as larger local building blocks in a lattice 
theory greatly enlarges the local space to be diagonalized in the beginning [26, 27] and 
enforces a restricted choice of states and/or the use of simpler impurity solvers like e.g. a 
simplified finite U- version of NCA (SNCA), especially if a self consistent determination 
of one-particle states and matrix elements is aimed at as part of the solution of a lattice 
problem [28] . (3) The numerical methods mentioned do not work equally well in different 
energy regimes. 

NRG, for example, often leads to a poor description at high excitation energies. 
Depending on tunable parameters in lattice models corresponding to temperature, 
pressure, chemical composition, doping, etc. a wealth of phases and corresponding 
transitions is found; even a transmutation of underlying pictures or concepts may occur, 
as e.g. from magnetism of (nearly) stable local magnetic moments coupled via short 
ranged exchange interactions to itinerant forms of magnetism to be described by an 
effective Stoner-theory for bands of itinerant quasiparticles [21, 29]. In these cases 
a good description of a large regime of excitation energies is desired, which requires 
impurity solvers, which work equally well in a broad range of energies and allow for a 
controlled approach to known limiting cases. Facing these difficulties for a more complete 
understanding and description of lattice systems, further developments of methods with 
analytical background seem necessary. 

In this paper we will present a version of direct perturbation theory for an 
(effective) impurity, which combines elements of some existing improvements of NCA 
and goes beyond them in some respects. An essential ingredient of this new stage 
of approximation are processes with crossing band electron lines, as depicted along 
a (imaginary) time axis, and hence the name crossing approximation one (CA1). In 
the following section 2 we will give a very short account of direct perturbation theory, 
describe CA1 and comment on its relation to NCA and its extensions as accounted for 
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in the literature so far. Section 3 contains some calculated spectra for the Anderson 
impurity model and a comparison to former approximations and, most important, to 
results obtained with the NRG. As it will become apparent, CA1 turns out to be a rather 
good impurity solver in the whole range of energies and can be applied to calculate 
excitation spectra, ionic propagators and susceptibilities. A conceptual bridge between 
impurity and lattice theories is outlined in section 4, which allows to recover the results 
of lattice theories like XNCA [20] and DMFT [22, 30] in a general fashion. 

As an example for the usefulness of CA1 in this context a calculation for the 
Anderson lattice is presented. The concluding section 5 contains remarks about an 
application of direct perturbation theory to susceptibilities and magnetic phases of 
lattice models and about improvements regarding nonlocal correlations as well as a 
new impurity solver, i.e. a CA2-project. 

2. Description of CA1 

2.1. Introductory remarks on direct perturbation theory 

A typical setup for the application of direct perturbation theory uses a Hamiltonian 
H = H £ + H 0c + V with the following parts: H £ = H^{{f m(7) contains the 

dynamics of interacting electrons in local one-particle states with quantum numbers 
m and a (pseudo spin) and is expressed via corresponding annihilation (creation) 
operators fmJ- H 0c = ^o c ({%j, c ka\) describes a reservoir of noninteracting electrons 
in Bloch states (a band index is suppressed here), and V = V{{f m(7) f+ a , c^, c^}) is a 
hybridization or transfer between local and band states, which is likewise expressed via 
elementary one-particle processes. H £ acts on a local Fock-space of finite dimension, 
typically one or a few valence shells or orbitals, and can in principle be diagonalized. 
A basis of eigenstates | n ,M ) ("ionic states") is denominated by a local particle 
number n and a set of many-body quantum numbers M specifying angular momenta 
or crystal field levels. With the operators X n / M t inoM =| n f M f )(n M | and corresponding 
n -particle energies E UoM the local Hamiltonian reads 

Hq£ E noM X no M,n M (1) 

where only projectors onto the eigenstates appear. The terms "local" or "ionic" do 
not necessarily imply one single atom. The formalism equally well applies to local 
subsystems of molecular type or to local clusters. A transcription of V to local many 
body states involves via 

fiS = E E a ^ " 1M '' n oM)X { ^_ 1M , noM (2) 

the set of ionic transfer operators X n /M f ,n M with n' Q = n ± 1. Using V as the 
perturbation, processes of direct perturbation theory are constructed from elementary 
absorption or emission events of band electrons from a local shell state at fixed 
(imaginary) times with amplitudes given by the coefficients a ma (n f M f , n$M) in (2) [3]. 
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Insofar it can be expected that e.g. the partition function can be cast into a form, where 
the contribution of each particular ionic state becomes apparent: 

Z = Tre-f> H e-f"Tr e Tr e {z - H)- 1 = Z 0c ^ <f ^ e^'P^z) 

= Z 0c ^ / duje~ Puj g noM (u). (3) 

Z 0c = Tr c e~ f3H ° c is the partition function for the band part alone and — 
— \lm PnoM^ + iS) is the spectral intensity of the ionic state |n M), which evolves 
with the propagator P no M(z)- A straightforward concept of irreducibility with respect 
to intermediate ionic states allows for the introduction of irreducible ionic selfenergies, 
with analytical properties as usual, 

Puom(z) = (z- E no M - S noM (^)) _1 , (4) 

and a corresponding perturbation expansion. The processes contributing to these 
selfenergies S no M(^) will in the following be constructed from skeleton diagrams, so 
that the ionic propagators P no M(z) are to be determined selfconsistently from a set of 
coupled integral equations. 

Representations for general Greensfunctions, which in correspondence with the 
partition function (3) are expressed as convolutions of ionic propagators, can also be 
derived along the lines sketched above. We consider in particular the local one-particle 
Greensfunction, 

F ma {r) = -(T(/ mCT (r)/+J) 

= -s~2Yl a ma (n -lM[,n M 1 )a* m(7 (no-lM^noM 2 ) 

(T(X no _ 1M ^ noMl (r )Xn M 2 fi -lM! i )), (5) 

the Fourier-coefficients of which at Matsubara frequencies cu n = ( 2n + 1 ) 7F ( n g Z) give, 
after analytical continuation iuo n —> u + iS to the upper border of the real frequency 
axis, the local one particle excitation spectrum Q ma {uS) = —^ImF ma (uj + iS). The 
complete setup of nonstandard direct perturbation theory is well documented, including 
the diagrammatic rules for the processes to be discussed in the following [3, 8, 31]. 



2.2. General features of crossing and non crossing approximations for the SI AM 

First comprehensive studies of direct perturbation theory for the Anderson impurity 
model (SIAM) 

h = E ( e * flf* + f ^) +E E { V AJ° + h - c ) 

concentrated on the limit U — oo of infinite local Coulomb repulsion and were termed 
NCA [8, 9, 32]. They were based on the leading skeletons of order V 2 to the ionic 
self energies T,q(z) and T,i a (z) and furnished a qualitatively correct picture e.g. for 
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the temperature dependent formation of the Abrikosov-Suhl resonance (ASR), the 
most prominent many-body signature of the Kondo-effect [33] in the local one-particle 
spectrum of the model. The particular aspect of degeneracy v of the local level, i.e. 
v = 2 for the two possible z— components of spin in the original SIAM but higher v 
as in Ce-compounds with v — 6 becoming possible through orbital degeneracy, drew 
much attention: In the limit v —> oc, using a proper scaling V —> the NCA-results 
become increasingly valid [2, 13, 34], and for v —> 1 on the other hand a trivially solvable 
resonant level model emerges. Whereas the limit of large v gave reason for classification 
schemes of diagrams in orders of ^, the limit v — 1 arose hopes of reconstructing an 
exact solution of a simple model by direct perturbation theory thus completely clarifying 
the systematics of diagrams for all cases of v. 

This hope was not fulfilled up to now, although Keiter presented an exact solution 
of SIAM for the zero-bandwidth limit unravelling the full diagrammatics for this simpler 
case [35]. It became clear then, that progress with the direct perturbation approach had 
to be worked out stepwise by including more important classes of skeleton diagrams into 
the calculations. 

In the following we discuss the systematics of these approximations for the SIAM 
by concentrating on the vertices, which allows for writing down several quantities in a 
compact and rigorous form, i.e. the ionic selfenergies 

dxD a (x)f(x)A 0M (z,x)P la (z + x), 

a J 

SiaO) = J dxf(x) D a (-x)A 0M (z + x, -x)P (z + x) (7) 

+ D_ (J (x)A 2 ,i (J (z + x, -x)P 2 (z + x) , 

£2(2) =^2 / dxD -cT(-x)f(x)A 2M (z,x)P la (z + x), 

and the local one-particle Greensfunction (5), which in the special case of the SIAM 
contains only two contributions: 

F a (z) = ±-J ^e-W [Ao,i.(^, z)P (z')P la (z + z') (8) 
c 

+ A 2| i_ <7 (^, -z)P 2 (z')P 1 _ a (z - z') 

Here we have introduced the hybridization intensity D a (e) = "}2 k S(e — e^) = 

V ( e ) 2 Qoc(t) and the (perturbed) local partition function Zn — 7^-. Expressions for 
higher Greens functions take an analogous form; we only add here a formula for the 
dynamical magnetic susceptibility (leaving out prefactors (^gfjip) ), 

1 P dz' 

Xma 9 {z) = - Y J^-e-^J2A a , (7 (z\z)P 1(7 (z , )P 1 Az + z'), (9) 
c a 



CONTENTS 



8 




Figure 1: Diagrammatic representation of ionic selfenergies, local one-particle 
Greenfunction and magnetic susceptibility for a single-impurity Anderson model (SIAM) 
in direct perturbation theory. Physical processes are arranged vertically along an 
imaginary time axis (broken line) which bears an energy variable z after Laplace- 
transformation. Presence of an electron in the local shell is indicated via a wiggly line on 
this time axis. Excitations of band electrons (straight lines) take place at hybridization 
vertices (dots on the time axis). Due to time-rotational invariance all vertex corrections 
in these diagrams can be collected at one of the vertices, which is drawn as a triangle. 

which involves a separate kind of vertex A^. Eqs. (7) to (9) are graphically represented 
in figure 1; observe identities like Aq^^ — z',z') = Ai a ^(z,z f ) for setting up the 
equations from there. 

In the SNCA, which can be viewed as the simplest nontrivial approximation for all 
values of £7, the vertex functions are all taken without any vertex corrections: 

SNCA: A 0M = h 2M = A^ = 1 . (10) 

The original NCA constitutes the U —> oc-limit hereof and is obtained by ignoring 
the doubly occupied state, i.e. by setting P 2 = 0. 

The NCA, and in consequence the SNCA for general values of {/, can only 
furnish qualitative insight into the dynamics of the SIAM, since it is plagued by 
shortcomings. These are revealed in the following ways: (1) Comparison with the 
resonant-level limit v —> 1 ends in an insufficient fit to the virtual scattering resonance 
in the excitation spectrum, in particular when the intermediate valence regime is 
approached [8]. (2) Form and position of the ASR are not in accord with the 
Friedel sum rule, most important for v — 2, and correspondingly the local self energy 
Ya g {z) = T l(7 (z) + iAA = z — 6£ + iAA — F a (z)~ 1 does not comply with local Fermi-liquid 
properties [36] (A^ = ttV (0) 2 Qo c (0) is the Anderson width and = E\ G — E the local 
one-particle level). (3) threshold exponents, as taken from the ionic propagators (see 
below) with values a = | and a\ a = | in the NCA, do not agree with the values 
known from the X-ray-absorption problem [17]. (4) An exact analytical solution of the 




Figure 2: Vertex structure of the "full NCA" (FNCA). The hierarchy of vertex 
corrections is generated by two coupled integral equations. The bare vertices ( first 
terms) are subsequently crossed by one more band excitation, which ends below and 
above in a full vertex, respectively. The diagrams are taken as skeletons, i.e. the local 
lines are dressed with the full ionic propagators. 

NC A- version of Eqs. (7) for zero temperature and a flat conduction band density of 
states symmetric around the Fermi energy reveals spurious features near the ASR [15], 
namely a sharp spike showing up at the Fermi level below a (" pathology-" ) temperature 
T p , being still lower than T K in the Kondo regime. 

Point (3) deserves some further comments, because it hints to the particular singular 
structure of the ionic propagators Pn M{z) ) which causes difficulties in the numerical 
solution of the system (7) of integral equations and also in subsequent procedures like 
(8), (9) involving convolutions of several of the P no M- As explained e.g. in [17], these 
propagators develop a common threshold at an energy uu = E g < for zero temperature, 
due to a slow algebraic decay P no M(t) ~ e~^ Egt /t aM in the time domain, i.e. 

2 2 

< Tin < 1 being the occupation of the local level in the Kondo-regime < < 
— €£ < U. At T — 0, E g is the lower endpoint of a branch cut in the functions P no M(z) 
along the real axis z = uj > E g \ it is this particular divergent structure - for < 1 and 
v — 2 one has 1 — a G > \ and < 1 — ol\ g < | - which needs care and makes numerical 
calculations to higher orders much more time consuming than NCA or SNCA, due to 
multiple convolutions of these singular structures. 

2.3. Generalizations of SNCA: PNC A, ENCA and FNCA 

The first useful generalization of NCA to the SIAM with general values of the Coulomb 
repulsion U was proposed and investigated in 1989 [37]. It was called "full NCA" 
(FNCA) and is visualized diagrammatically in figure 2. One recognizes a particular 
subsystem of integral equations, which serves to generate a class of vertex corrections 
(again as skeletons) extending to infinite order. This particular choice was motivated 
by an attempt to include as many as possible exchange counterparts to those processes, 
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Figure 3: Sequences of two elementary excitation processes, which together 
constitute the lowest order exchange coupling vertex remaining after a Schrieffer- Wolff 
transformation of SIAM to the s-d-model. 

which already contribute to the ionic propagators in SNCA, see appendix in [37]. As 
Keiter repeatedly has pointed out [38], the balance between processes which transform 
into each other by a reversal of partial time orderings, as shown in figure 3, is necessary 
to obey the Pauli-principle and to comply with universality in the Kondo limit, where 
in accord with the Schrieffer- Wolff transformation from SIAM to the s-d-exchange 
model [39] the characteristic energy scale k B T K is expressed via an effective exchange 
coupling constant J = — ^ttj \ figure 3 just visualizes the two contributions to this I 



The system of five integral equations according to figures 1 and 2 was solved for finite 
U in [37], and the results for the ionic propagators and the corresponding excitation 
spectra were compared to some simpler calculation schemes. Whereas pronounced 
discrepancies to SNCA showed up, e.g. regarding the important energy scales, the so 
called "enhanced NCA" (ENCA) already captured important improvements. 

In ENCA all vertices on the right hand side of the two equations in figure 2 are 
taken as bare ones. Then, only the leading contributions to the infinite series of vertex 
corrections contained in FNCA are included; among the latter are running n-particle 
cascades between initial and final state during the excitation by the external electron 
(iterate the vertex in the middle of the last diagram) as well as long-time memory effects 
between initial and final states through chains of internal excitations, before or after the 
external excitation occurs (iterate the respective vertices on top and at the bottom of 
the diagram). Since the ENCA has proven as a good compromise between accuracy and 
the calculational effort to be invested in an impurity solver for lattice problems (see also 
section 4), we cite the explicit expressions for the vertex corrections, which have to be 
solved together with the system (7) of self-energy equations: 



Calculations of the local one-particle spectrum in [37] were then based on the ENCA 
and led to an improved many body scale and a better understanding of the many body 
dynamics of SIAM, in particular at finite values of U. 



[11, 40]. 




(z, z') = j deD_ a (e)f(e)P 1 _ a (z + e)P 2 (z + z' + e), 
(z, z') = [ deD. a (e)(l - f(e))P la (z - e)P (z - z' - e). (12) 



^^2,1-0- 
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Figure 4: Part (a) shows a fully crossing vertex correction of order 0(V 4 ). The two 
processes shown in parts (a) and (b) are of the same order regarding an expansion in the 
degeneracy v of the singly occupied local state. With the self consistent replacements 
shown in part (c) and with fully dressed local lines they constitute the "post-NCA" 
(PNCA), a theory for U = oo, in which the doubly occupied local state is projected out. 



Whereas ENCA takes into account the vertex corrections up to order Oiy 2 \ and 
FNCA in addition certain classes up to infinite order, both do not include the fully 
crossing diagram of order 0(V 4 ) shown in figure 4(a). This vertex correction is the 
lowest non vanishing one in the U = oo -theory and was frequently used to discriminate 
"crossing" and "non crossing" approximations. 

In order to investigate the role of such fully crossing diagrams the SIAM at infinite 
U was investigated in 1994 with help of a "post-NCA" (PNCA) [41]. This approximation 
scheme was set up along the lines of a —expansion and collected all vertex corrections 
up to O (^2), i.e. all contributions to the ionic self energies up to this order. Therefore, 
also the vertex correction shown in figure 4(b) was taken into account, which has two 

more powers of V compared with figure 4(a), but due to [~%) = ^~ ^ s °f the same 

order ^ thanks to the closed ring with spin-summation over a f between vertices 2, 3, 6, 
and 7. Actually, and in close analogy to the FNCA, the bare vertices in figure 4(a) and 
(b) were all replaced by full ones, as indicated in figure 4(c), and the coupled system of 
vertex corrections (now including all orders) was solved, again self-consist ently together 
with the system (7) of ionic self energies. Convergence could be reached on not too large 
time scales by use of parallel computing. Progress over the original NCA turned out 
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to be essential: Apart from a corrected many-body scale k B T K , the local Fermi-liquid 
properties improved considerably, the position of the ASR near the Fermi level agreed 
much better with the one implied by Friedels sum rule, and also the threshold exponents 
a and a\ were shifted towards the values of (11), although agreement with these values 
or with a variant according to [42] was not conclusive. 

Due to the considerable numerical effort, regarding the multiple overlapping 
integrations over functions with rich structure, an extension of PNCA to finite values of 
U seemed not possible in 1994, since many more diagrams involving the doubly occupied 
state would have to be added. 

2.4- Other approximation schemes in the literature 

Before the new approximations CA1 and a CA2-project (in section 5) will be explained, 
we shortly comment on two approximation schemes, which have been proposed and 
investigated over the last ten years. In the so called "symmetrized finite-t/ NCA" 
(SUNCA) special emphasis is laid on the chains of scattering events [43] mentioned 
above in connection with the FNCA [37]. This scheme is conserving (^-derivable in the 
sense of Kadanoff and Baym) like all other approximations mentioned in this section; 
moreover it can be characterized as involving just a subclass of the FNCA-diagrams. 
Although the relevant papers are written with help of the slave-boson formalism, the 
formulation is fully equivalent to direct perturbation theory as pointed out above. The 
evaluation of the local one particle spectrum is based on a full infinite subclass of vertex 
corrections and thus goes beyond the ENCA-calculations. These vertex corrections 
with long scattering chains are easily visualized with help of the FNCA-diagrams of 
figure 2: Iterate the vertex equations with respect to the upmost vertex only. The 
results underline the progress reached with ENC A and FNCA [37] . 

Whereas SUNCA is applicable to the finite- U case and can be placed into a scheme 
of repeated vertex corrections with single line crossings (a more general version being 
FNCA), the "conserving T-matrix approximation" (CTMA) [44] again is restricted to 
infinite U and stresses the importance of chains of scattering events for band electrons 
off the local shell over the whole duration of the external excitation process. These are 
argued to contain those significant contributions, which are known to lead to the correct 
singular threshold behaviour of X-ray absorption spectra as predicted by Mahan [45] 
and calculated by Nozieres et al. [16]. 

Correspondingly, essentially exact threshold exponents are expected from the 
CTMA. This approximation can be characterized with reference to the fully crossing 
diagram of figure 4(a): The middle part between vertices 2 and 3 becomes the lowest 
contribution to a T-matrix, which is fully determined by the implicit equation shown in 
figure 5(a). It generates the sequence of vertex corrections with scattering chains shown 
in figure 5(b). Observe that only the first of these is contained in PNCA. 

In spite of a superficial resemblance already the second contribution is different from 
figure 4(b), which is more easily recognized by counting the number of independent spin- 



CONTENTS 



13 




Figure 5: Vertex corrections summed in the "conserving T-matrix approximation" 
(CTMA) for the SIAM at U = oc as viewed from direct perturbation theory. Part 
(a) T-matrix which, when substituted into the diagram of figure 4(a), generates the 
sequence of vertex corrections shown in part (b). 



summations. Indeed, CTMA-results [46] point to considerably improved values of the 
threshold exponents; nevertheless, the description of the local Fermi-liquid formation, 
similar to PNC A, is still not fully satisfactory. Both of these approximations involve 
time-consuming numerical calculations; up to now a generalization to the even more 
demanding case of finite [/-values has not been reported. 

Other approximation schemes involving additional simplifying assumptions for the 
ionic propagators and vertex functions, be it either in a non-conserving [47] or conserving 
fashion [48], will not be considered here. Although they may be useful with respect 
to computational effort, they have only been justified for the case of large orbital 
degeneracy. 

2.5. CA1 approximation 

CA1 is designed to describe SIAM in the full range of values for the local Coulomb 
repulsion U with good accuracy and likewise for dynamical properties at general 
excitation energies uj. Being a straightforward collection of all vertex corrections up 
to order 0(V 4 ) (as skeletons) it is conserving and contains the leading contributions 
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Figure 6: CA1 collects all vertex corrections for general (finite and infinite) values of U 
up to order 0{V A )] in part (a) these are shown explicitly for one of the two vertices. 
The analogous construction for the other vertex is indicated in part (b) by dots. Again 
local lines are dressed, i.e. the diagrams are used as skeletons; the vertex points on the 
right hand side, however, are bare ones. 

from all of the approximations sketched above. More precisely, it can be defined by 
the set of vertex corrections shown in figure 6(a), plus the corresponding ones for the 
other vertex, indicated by points in figure 6(b). In successive order the diagrams may 
be characterized as follows: The first terms on the r.h.s. are the bare vertices and 
define SNCA. ENCA additionally contains the following vertex correction with a single 
crossing electron or hole line, respectively. The next three vertex corrections can be 
viewed as originating from the ENCA-diagram by dressing each of the vertices with 
a single crossing line successively, i.e. they represent the first iteration in the FNCA- 
scheme. The last diagram is the fully-crossing contribution not contained in the FNCA; 
it is the leading vertex correction in both approximations for infinite J7, PNC A, and 
CTMA. 

Whereas CA1 is explored in the following together with the other approximations 
mentioned, a CA2-project will be designed to add more vertex iterations like those 
included in FNCA and longer scattering chains like those of the CTMA, which can be 
incorporated via a T-matrix formalism. 
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3. Results from CA1 and comparison with other impurity solvers 

3.1. General remarks on the quality of impurity solvers 

The following criteria have frequently been applied to judge the quality of impurity 
solvers in connection with the Anderson impurity model: (1) Ionic propagators have to 
obey the correct threshold behaviour, in accord with the relevant work on orthogonality 
catastrophy and excitonic correlations in the case with spin-degeneracy [17, 45, 49, 50]. 
(2) The infrared divergencies of the perturbation series produce a characteristic low 
energy scale, usually referred to as the Kondo-temperature T K) which should faithfully 
be reproduced by the approximation. (3) The many-body resonance (Abrikosov-Suhl- 
resonance, ASR) forming at temperatures of order of the Kondo-scale T K and lower is 
pinned at a position near the Fermi level, which is determined by Friedels sum rule. (4) 
The shape of the ASR has to comply with a form of the selfenergy, which guarantees 
local Fermi-liquid properties. 

In the following, the four semianalytical impurity solvers introduced hitherto for 
the SIAM with general, in particular finite values of U, i.e. SNCA, ENCA, FNCA, and 
SUNCA will be compared to the new CA1. Special emphasis will be laid on the above 
four criteria. As a reference, also NRG-calculations for the SIAM are presented, which in 
the low-energy regime should provide a reliable bias. They have been produced with help 
of the very effective numerical procedures presented in [51]. All other calculations have 
been performed with a software package written for the solution of a number of impurity- 
and lattice-models, in which various impurity solvers can be combined with different 
methods for the lattice aspects. Since it is based on adaptive strategies for an all- 
purpose use, no particular provisions have been taken to optimize numerical strategies 
for the SIAM in the deep Kondo limit. Nevertheless, the program package seems to 
work very reliably, although numerical convergence problems and approximation errors 
become visible for certain extreme choices of model parameters. In particular, the 
folding of several nearly singular factors in an integrand like that of (8) at very low 
temperatures, needs a thorough analytical preparation and consumes much numerical 
effort, and likewise the iteration of vertex parts depending on two energy- variables for 
SUNCA, FNCA, and CA1. 

3.2. Ionic threshold behavior 

The shortcomings of SNCA have already been mentioned in the foregoing section. It 
proves worthwhile, however, to check how the known values of threshold exponents for 
this approximation are recovered in the calculations; this gives valuable hints as to how 
results for the other approximations should be interpreted and generally, how reliably 
the algorithms work. It is interesting to note here, that the information given about 
the NCA-threshold exponents (i.e. the case U = oo) in connection with (11) is not 
complete when regarding SNCA at finite U. The treatment of [15] based on an ansatz 
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Figure 7: Double-logarithmic plot of ionic spectra, centered at the threshold, for a 
SIAM in SNCA (a) in the symmetric case with e t = -1.0, U = 2.0, f3 = oo, (b) in 
an asymmetric case with q = — 1.0, U — 3.0, (3 — oo, both for an Anderson width 
= ttV 2 gf} (0) = 0.3 and a 3d — sc band density of states Qw(uo), centered at 
uj = ji = 0; the bandwidth is 6 units. The threshold exponents can be read off as 
one plus the slope of the asymptotic tangents drawn in the Figures. 



for ionic self energies is easily generalized: 

S M (^ + iS) ~ Eg — Em — - Eg) 1 '*™ P noM (u + iS) « -iA M (u - E g ) aM -\l3) 

It furnishes the same asymptotic region near u = E g as in the case U — oo, characterized 
by a = | , alcr = | for (spin-) degeneracy two, for the whole unsymmetric regime 
(c£ = Eio- — Eq < —A and) 2e£ + U > 0. Asymptotically, the propagator P 2 does not 
contribute here. For U approaching the value — 2q from above, however, P 2 becomes 
equal to P , which leads to 

a i = 2^' ao = a2 = ^; (2e/ + t/ = 0), (14) 

being the degeneracy of ionic state |la). For the model with spin-degeneracy only 
this means — a\ — a 2 — |. These values coincide with the presumably exact ones 
taken from (11). Naturally this does not imply that SNCA becomes correct for the 
symmetric SIAM. As we will see below, e.g. the shape of the ASR still reveals serious 
shortcomings. The transition from U — oo to U = — 2q happens in a gradual way: 
The former asymptotic regime around uu = E g shrinks to zero, and a different regime 
takes over, which originally was situated at higher values of \uo — E g \ and developed the 
exponents of (14). 

In figure 7 and 8 we show the results of a rather precise SNCA-calculation of a 
symmetric SIAM in the deep Kondo regime, parts (a), and an asymmetric one, parts 
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(b). The impurity states locally hybridize with a tight-binding simple cubic conduction 
band in three dimensions of width 6, Fermi level and band center lie at energy uj = 0, 
and van-Hove-singularities at uj = ±1. Shown in figure 7 are the spectra of the relevant 
ionic propagators P and P\ G for T = on a doubly logarithmic scale with origin at 
the corresponding threshold energies E g . Numerical resolution is somewhat below 10~ 4 , 
acting as an effective temperature cutoff. The asymptotic regime is entered only one 
order of magnitude higher, at about uj — E g « 10~ 3 . From the slope of the tangents 
drawn one reads off the exponents a ~ 0.53, a\ ~ 0.47 in the symmetric case, and 
a = 0.38, ai = 0.76 in the asymmetric one. Even if by a proper extrapolation, using 
several low values of the temperature, these numbers can be brought closer to the exactly 
known ones, i.e. (0.5, 0.5) and (0.33, 0.67), an uncertainty of about 10 percent remains. 
This is enough, however, to identify and discriminate the two different situations. With 
lower numerical accuracy or at higher temperatures T > 10" 4 « TV/200 the asymptotic 
regime will be hard to attain. 

Figure 8(a) demonstrates that in spite of accurate threshold exponents the shape of 
the ASR at uj = in the 1-particle-spectrum for low T comes out as a quite unphysical 
spike. Part (b) of this figure shows how the ASR is deformed, when U is raised; the 
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position of the ASR is clearly temperature dependent, and its flank at uo = develops 
a pathological steepness. This reminds of the pathological structure found at infinite 
U with a flat conduction band density of states [15]. Since for our calculations a three 
dimensional simple-cubic tight-binding bandstructure was used, van-Hove singularities 
(i.e. kinks at uj = ±1) leave their traces in the spectrum; their visibility also constitutes 
a test for the numerical procedures used. 




Figure 9: Spectral density of the empty ionic state M = for a SI AM, calculated 
with the five semianalytical approximations discussed in the text, parameter values as 
in figure 7, (3 = 1000. The shifted thresholds allow for fits of different quality to the 
Kondo temperature T K . 



A good qualitative insight into the relation between the five semianalytical impurity 
solvers under consideration can be obtained from figure 9, which shows the spectrum of 
= P2 f° r the symmetric SI AM discussed before at temperature T = 10 -3 . Whereas 
part (a) gives an overall view with the threshold to the left and a broad one-particle 
resonance to the right, corresponding to a distribution of contributing frequencies around 
uj = —Ei a —AEi a (E is set to zero in all calculations), part (b) with a much finer energy- 
resolution points to the discrepancies between the different approximations visible in 
the low-energy regime. One recognizes threshold-peaks at different values of E g) in 
increasing order for FNCA, CA1, SUNCA, ENCA, and SNCA. Differences in E g directly 
reflect the ability of the approximations to reproduce the Kondo-scale, which in this 
regime can be expressed as [52] : 

[k B ]T K = aVlexp\-j} , 1 = . (15) 

L 1 J eAep + u) 



Choosing a — ^ for U < bandwidth W [37] we obtain T/T K = 0.04, i.e. the spectra 
represent the temperature range well below T K) even with slightly different choices of the 
coefficient a. If one accepts for the moment, that the FNCA with lowest E g furnishes 
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Figure 10: Double-logarithmic plot of ionic spectra, centered at the threshold, for a 
SIAM in ENCA, parameter values as in figure 7. 



the closest approximation to the real T K) as e.g. is implied by the NRG-calculation 
(see below) then the following conclusion can be drawn: ENCA, SUNCA, and CA1 all 
improve considerably on the SNCA. The leading vertex correction already included in 
the ENCA contributes the essential part to this effect, whereas the additional terms 
further taken into account in SUNCA and CA1, respectively, have a relatively smaller 
impact. This agrees with the original investigation of ENCA and FNCA [37], where 
it was shown, that the ENCA already captures the right exponential behaviour of Tk 
for the SIAM at finite U ) whereas the inclusion of further vertex corrections then only 
improves on the prefactor in this scale. 

As will be shown below, the good estimate of T K furnished by the FNCA does 
not imply that FNCA behaves well in all other respects, e.g. concerning the four 
points mentioned in the beginning. What is apparent, however, is the pronounced 
and qualitatively similar threshold behaviour visible in all of the five approximation 
schemes applied to the ionic spectra. 

Since ENCA, SUNCA, and FNCA all reduce to the NCA in the limit > e £ 
fixed, U —> oo, it is to be expected that in the asymmetric case, 2q + U > an 
ultimate asymptotic regime very near to cu = E g exists, where the ionic spectra are 
ruled by the NCA-threshold exponents. This must not necessarily be true for the 
CA1 with its fully crossing vertex correction, which does not vanish in this limit. It 
can nevertheless be anticipated, that as a precursor a regime with "better" threshold 
exponents at somewhat higher values of \uo — E g \ occurs also for ENCA, SUNCA, and 
FNCA. In figure 10(a) a corresponding evaluation of ENCA is shown for a temperature 
which again is far below T K and also below the numerical resolution of about 10 -3 or 
somewhat less. The exponents for the symmetric case in figure 10(a), as read off in 
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Figure 11: Double-logarithmic plot of ionic spectra, centered at the threshold, for a 
SIAM in SUNCA (part (a)), FNCA (part (b)) and CA1 (part (c)), all for the asymmetric 
case with parameter values as in figure 7. 



the range 10 -3 < u < 10 -2 are clearly near the value 0.5, whereas in the asymmetric 
case of figure 10(b) only a « 0.5 is really conclusive; a value ol\ ~ 0.5 can be justified 
only if a tangent is drawn in the reduced range between uj ~ 10 -3 and uj ~ 0.5 • 10 -2 
before the steeper decrease sets in. Comparing figures 7(b) and 10(b) this seems to 
be a reasonable procedure. It must be remembered here that the better approximation 
cannot be evaluated with the same numerical accuracy, at least not using the program 
package in its present form. Figures 11 (a), (b), and (c) show the threshold behaviour 
for the asymmetric model (again with e\ = — 1.0, U = 3.0 and T <C T K ) obtained with 
SUNCA, FNCA, and CA1. Similar conclusions as for the ENCA can be drawn here: In 
all cases the threshold exponents in the accessible asymptotic regime come out near the 
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Figure 12: One-particle excitation spectrum of a symmetric SI AM at f3 — 100, other 
parameters as in figure 7(a), in the five semianalytical approximations discussed in the 
text and, additionally, calculated with the numerical renormalization group (NRG). 
Part (a) reveals shortcomings of the NRG-method at large excitation energies, whereas 
in the low-energy region of part (b) the NRG-curve can be used as a reference for the 
other approximations. 
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Figure 13: Analogue of figure 12 for the asymmetric case with f3 = 150, other parameters 
as in figure 7b. 

value 0.5. At least for the CA1 this gives reason to hope for an essential improvement 
of the true asymptotics of ionic spectra over the SNCA. 
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Figure 14: Like figure 12, but with (3 = 1000, i.e. at T < T K . 



5.5. One-particle excitation spectra 

For a comparison of the one-particle excitation spectra obtained with the five 
approximation schemes also a calculation with the NRG is taken into account. This 
should give an impression of the exact result, at least in a low-energy regime near 
uj = fi = 0. 

In figure 12 we present the results for the symmetric SI AM and in figure 13 for 
the asymmetric one, at (3 = l/fc^T = 100 and (3 = 150 respectively, implying T being 
roughly half T K) other parameters as before. First looking at parts (a) of these figures, 
the following fact seems remarkable: Peaks in the NRG-spectra are considerably broader 
compared with the other cases and the features due to the van-Hove singularities in the 
band-DOS are smeared out much more. At least part of this, in particular at higher 
excitation energies, should be due to numerical procedures used in the NRG-calculation: 
It is based on a discrete set of energy-eigenvalues, which is considerably spaced near 
the bandedges and becomes logarithmically denser for smaller excitation energy uj; 
interpolation then has a smoothening effect. Whereas near the original resonances at 
uj < 6£ and lu > 6£ + U the five semianalytical approximations are certainly closer to the 
truth than the NRG-curves, the situation is not completely clear in the energy region 
around the ASR, although here the NRG is most trustworthy. 

The peak value close to uu = at T = is given via Friedels sum rule to be 
Qfa{uS) = 1/ttAa ~ 1.06; this is faithfully reproduced by the NRG, whereas the width 
of the ASR might already be somewhat exaggerated by the NRG. In effect, however, 
we take the NRG- ASR as our measure of quality for the other approximations with 
respect to the low-energy regime. Regarding the full range of excitation energies, on the 
other hand, any of the other approximations (except SNCA) might be more appropriate, 
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depending on the purpose of the calculation. With the halfwidth of the ASR taken as 
measure for the many-body scale T K) certainly the FNCA, with its vertex-corrections 
systematically iterated through all orders, compares most favourably with the NRG. 
On the other hand, the FNCA-spectrum clearly exaggerates the height and thus the 
total weight of the ASR: The limit of 1/ttAa for the peak-height becomes violated even 
stronger for lower temperatures. Insofar, CA1 seems to represent a good compromise 
and even ENCA does not work too bad. 

The lesson to be learned from these results is that improved semianalytical impurity 
solvers of this type should incorporate skeleton diagrams of two types in a well-balanced 
way: Classes of iterated vertex corrections have to be accompanied by chains of iterated 
particle-scattering events, being related to each other as exchange-partners [37, 44, 46]. 
In light of the discussion in section 2, CA1 serves as a further step in this direction. 
Regarding the regime at high excitation energies uj ~ q + U the comparison of 
figures 12(a) and 13(a) reveals a trend, which for even larger values of U becomes 
more and more pronounced and which apparently is not well captured by the NRG, at 
least with its present numerical performance: The resonance due to double occupancy 
of the local shell becomes sharper with increasing £7, in particular when e^ + U reaches 
the order of the upper band edge uj « \W = 3. Beyond this value, the peak keeps 
its weight but rapidly looses its width and finally vanishes as a single spectral line out 
of the accessible region. This is faithfully reproduced by any of the five semianalytical 
impurity solvers under consideration. 

Figure 14 gives an impression about qualities and failures of the five approximations 
as applied to the full calculation of one-particle spectra at very low temperatures: 
= -1.0, U = 2.0 and (3 = 1000, i.e. T/T K = 0.04 have been chosen here. Whereas 
the ionic spectra in all five cases come out rather reliably with the procedures used in our 
program package, the subsequent folding of ionic propagators and defect propagators 
(see e.g. [14]) can produce spurious results near uj = 0. With very sharp thresholds in all 
quantities at low T slight displacements of the maxima (as a consequence of numerical 
procedures and rounding errors) can have a large effect on the integrals containing 
several of these quantities. Although figure 14, too, supports the conclusions drawn 
before, the SUNCA-curve and to a somewhat lesser degree the FNCA-curve, show a 
spurious double-peak structure near uj = 0, supposedly due to such threshold-shifts. 
In addition, the FNCA-curve should not be taken too seriously very close to uj = 0, 
although its shape is in accord with the numerically more precise SNCA-calculation in 
figure 8(a). FNCA overestimates the ASR-peak height strongly, whereas CA1, in spite 
of a too high peak value, rather favourably compares with the NRG-curve. 




Figure 15: Imaginary part of the selfenergy of local electrons (absolute value) for a 
SIAM, calculated within CA1, as test for local Fermi-liquid properties. Temperatures 
as shown, other parameters as in figure 7. 



3.4- Fermi-liquid properties 

As a final test for our new CA1 the local Fermi-liquid properties are inspected in 
figure 15. For this purpose the imaginary part of the selfenergy 

-Im V fa («, + iS) = -— f \% a ^ ™ (16) 

is shown for a few temperatures near and well below T K . The formation of a minimum 
at uu = obviously takes place. In the asymmetric case of figure 15(b) a displacement of 
this minimum away from the Fermi level with growing temperature is recognized, similar 
but somewhat weaker than has been reported before for (S)NCA and ENCA [37], as well 
as for the PNCA, the latter being a U = oc-theory with crossing contributions to vertex 
corrections in high orders [41]. The CAl-result for the value of -Im T l f a (u + iS) at its 
minimum falls short of the exact limiting value (^^(a;)) -1 = = 0.3 for the lowest 
temperatures. This is in accord with the too high ASR-value seen e.g. in figure 14. In a 
quantitative sense, however, this result improves on ENCA and, much more, on (S)NCA. 
Obviously, the minimum can well be fitted by a parabola as long as the temperature 
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does not become so low, that numerical deficiencies near uo = become predominant. 
Also its position and height, as well as the quadratic coefficient may be compared to 
exact results for Fermi liquids [11, 58], e.g. as a guideline for corrective measures when 
using these impurity solvers in lattice-calculations. These conclusions are similar to 
those for the PNCA published before [41]. Furthermore, it is a remarkable fact how 
dramatic the scattering rate raises for increasing excitation energies. In the range of the 
ionic resonances u = and uu = + U it becomes high enough to completely prevent 
locally a band picture even for the c-electrons. This will become even more evident in 
the lattice calculations of the next section, where the effect occurs on each lattice site 
and thus affects the whole Bloch-states. 

4. Impurity solvers and lattice theories 

4-1- Cumulant perturbation theory of strongly correlated lattice models 

The importance of impurity solvers for approximate solutions of lattice problems came 
to light in theories like ATA [19], LNCA [21] and XNCA [20], which all used a picture 
of effective sites and relied on NCA as the best available implementation in those days. 
These three forms of effective site theories aimed at a solution of the Anderson lattice 
model in the context of the Heavy Fermion problem, and thus were based on a particular 
local shell structure with well localized /-states and extended c-states, with only the 
former being subject to a local Coulomb repulsion. In most cases, transfer also was 
restricted to the c-states only, which together makes possible a reduction from the matrix 
formulation, envisaged in section 2 and shortly outlined below, to a scalar formalism. 

The three theories differed in the way, in which the dynamics on general lattice sites 
influenced the representative effective site considered: In ATA only the coherent build 
up of quasiparticle bands from scattering by independent Anderson-impurities on the 
lattice sites was taken into account thus ignoring exhaustion problems and important 
renormalization effects. LNCA used self-consistently modified local excitations at the 
effective site; unlike XNCA, however, it introduced weight factors for nonlocal processes 
in order to approximately factorize the partition function into contributions from an 
unperturbed band and independent effective sites. XNCA finally established the form 
of self-consistency between effective site and surrounding medium which becomes exact 
in the limit of infinite spatial dimension and which nowadays is regarded as characteristic 
defining feature of the DMFT. 

In a local approach a lattice Hamiltonian H = Hq u + V y y is built up 

from local Hamiltonians H 0l/ on lattice sites R v , each of the type H £ considered 
in section 2, and nonlocal parts V v y, which contain one-particle terms like transfer 
or hybridization and possibly two-particle terms, i.e. nonlocal interactions between 
electrons on two different sites v and i/. V v y is usually expressed via the elementary 
creation and annihilation operators fltna which define the local Fock-spaces, whereas 
the diagonalization of H 0l/ involves the "ionic states" as described in section 2. In 
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principle there is complete freedom in defining what should be such a local subsystem, 
a single ionic shell, or an ion or a molecular complex, or even a cluster of ions or 
complexes. The need to diagonalize them, even when restricted to e.g. the low energy 
regime, may however set limitations to the size of what can be regarded as local. At 
the outcome will be a finite matrix in the space of local many body eigenstates 
\vtiqM). A one-particle transfer term contained in V v y will e.g. induce changes 
|z/n Mi) \vn -lM 2 ), \i/n' Q M[) -> \i/n' Q + lM£. 

A convenient way of keeping the formalism simple is to work generally with 
the original one particle quantum numbers (to, a) and to built up local matrix 
Greensfunctions of the type QL(z) = (Gm^mV'^)) f° r one-particle propagation, and 
corresponding higher ones; these matrices involve information about the composition of 
the local many-body eigenstates, e.g. via the coefficients of fractional occupancy. Thus 
an overall (N x A^)-matrix formalism for lattice processes is established with N being 
the number of local one-particle states taken into account, and local excitations of single 
electrons are translated via (2) into the dynamics of many-body eigenstates. Nonlocal 
interactions can be viewed as (simultaneous) two-particle transfers and be handled in 
an analogous fashion. 

The local approach treats the local dynamics exactly, as apparent in the many- 
body eigenstates. In the presence of strong and dominating local interaction matrix 
elements this seems a natural starting point for a lattice theory. However, it turns to 
a disadvantage when a perturbation expansion in terms of the nonlocal parts V v y of 
the Hamiltonian is to be set up. An attempt to factorize e.g. a general contribution 
to the partition function with the help of Wicks theorem stops at an intermediate 
level: In case of e.g. a pure one-particle nearest-neighbour transfer mechanism, time 
ordered expectation values on products of creation and annihilation operators belonging 
to sites, which knot together different propagation paths of single particles, remain as 
unfactorized parts of internally connected contributions B\. These also contain matrix 
elements t of V v y, symmetry factors r\, a sign (— 1) XA , site-summations and time- 
integrations, as well as an indicator function F assuming values or 1, which realizes 
site-exclusions between all the B\ occuring via n\ — 1). Schematically: 
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Figure 16: This contribution to the partition function is a product of three disconnected 
pieces, each of them containing single-particle loops. In two of the pieces loops are glued 
together at sites (nodes), which are marked with a double circle. These nodes give rise 
to cumulant vertices for local two-and three-particle interactions, respectively. 

In general, local correlations between all of the excitations caused at the nodes by 
intersite transfers remain. The emerging picture is that of disconnected sets of loops 
over the lattice, each set being internally glued together at certain sites with four or more 
intersite transfer legs, see figure 16 for a simple example. Moreover, while keeping the 
topological structure of such a graph, the position of sites involved cannot be summed 
freely over the lattice, thus preventing a convenient momentum space formulation. As 
a consequence, also a linked cluster theorem is not available for Greensfunctions of the 
lattice problem, since a partition function factor, the diagrams not linked to external 
sources, cannot be factored out: the partition function can also not be represented as 
an exponential of a sum of single connected graphs thus preventing a straightforward 
extensitivity property of the thermodynamic potential. 

A way out of this dilemma uses a representation of higher order local time ordered 
expectation values at nodes as a sum over products of successively smaller ones, which 
in total cancel out except for the original highest term; when appropriately grouped 
together they define a cumulant-expansion of local n-particle Greenfunctions, containing 
a set M n of n local destruction operators and a set M n of creation operators: 

G n {M n -M n ) = E (-i) Xp+x? U GC nS N S,$S)> 

allpartitions P of M n and P of M n into subsets (jSf^ q 

G\ (/ m ayD (r^), / m ( 2 y 2 )(7f ); f\ 2 y f\ iy (1) ,(r«')) = 

G 2 (l,2;2',l')-[G 1 (l;l')G 1 (2;2')-G 1 (l;2')Gi(2;l')] . (18) 

The right hand side of this expression can be viewed as containing two contributions 
playing each a different role in the expansion: The maximally decomposed terms, being 
products of only one-particle Greensfunctions at this site, together just furnish the 
result which would be obtained if Wicks theorem were valid for the local dynamics. 
The rest of the terms represent all possible local decompositions of the knot in the 
graph into products of independent knots of lower order with together the same number 
of intersite legs as the original expectation value; among them are possibly local one- 
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particle Greensfunctions but at least one of higher order, i.e. with more than two intersite 
legs. 

An evident approximation can be obtained as follows: Keeping only the maximally 
decomposed terms, all loops of the diagram, which were glued together originally, would 
be disconnected, the partition function would become an exponential of the sum of all 
different loops, and a free summation of all site positions over the lattice would leave only 
contributions differing by the number of links. Likewise, e.g. the lattice-one-particle- 
Greensfunction would contain, after dividing out the partition function, a sum over 
paths differing only by the number of links, which can be summed as a geometric series 
(in matrix space). 

The results for one-particle Greensfunction and thermodynamic potential F can be 
given explicitly as 

GT{z) = \G?\z)- l -t^ , 

F FT - F (0) = ~ Y. E^^f T) (^n)] - Tr[lna (0) (^n)]) , (19) 

with (2^ (z) the Greensfunction from the known solution of the isolated local subsystem. 
The above solutions may be viewed as a very general form of the well known Hubbard-I 
approximation [53] , to which they reduce when the local subsystem is a simple s-shell and 
V contains only nearest-neighbour one-particle transfers i.e. in the case of the Hubbard 
model. 

The rest of the terms in the cumulant-decomposition furnishes all possible 
combinations of one-particle scattering events off this site, serving as nodes with two 
intersite legs in simple transfer chains, and of higher order nodes with more intersite 
legs; together they realize a restricted way of glueing together the original connected set 
at this site. The contributions resulting from all of these local cumulant decompositions 
in the set can be re-interpreted in terms of a perturbation expansion with respect to an 
infinite set of local n-particle cumulant interactions, with n between 2 and infinity. 

The "non-interacting" starting point of this expansion is the aforementioned 
generalized Hubbard- I-theory free of such interactions ("Free Theory"). A perturbation 
expansion with the corresponding "free" propagators for the particles and the set 
of all local cumulant vertices as interaction terms along the conventional lines a la 
Feynman thus faithfully produces all contributions to partition function and lattice 
Greens functions and re-introduces the applicability of the linked cluster theorem and 
all benefits connected with it into this new form of the theory. 

The prize paid for the conceptual progress described above for the lattice problem, 
i.e. the applicability of conventional methods in perturbation theory for the lattice 
aspects, lies in the large number of vertices appearing and in their dynamical nature, 
i.e. their multiple time-dependencies. It should also be clear that the perturbation 
series obtained in this way are based, although looking conventional in a superficial 
way, on quite unconventional definitions of connectedness and irreducibility: Whereas 
the diagrams remaining i.e. for a one-particle-Greensfunction are linked to the external 
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sources and consist of one connected piece (we exclude anomalous terms here), being 
glued together via local cumulant vertices, they would in most cases fall apart into 
several unconnected pieces, only one of them bearing the two external links, when the 
cumulants are made explicit in terms of the original expectation values. 

The pieces without links remaining in this case as factors are partially due to the 
partition function in the denominator of the original expression and must consistently 
be included for a proper renormalization of e.g. excitation spectra. It is therefore not a 
trivial problem to define consistent approximations in cumulant perturbation theory. 

A straightforward evaluation of cumulant vertices even in low orders requires some 
effort but can be managed e.g. with methods of direct perturbation theory [54]. The 
somewhat lengthy expression of the local two-particle cumulant has been published 
before [55] and served as basis for calculations of one-particle properties of the Hubbard 
model along the lines of Hartree-like expansions [56] , also in a self-consistent fashion [54] . 
It should be clear and can in fact be proven [54], that these approximations are not 
able to capture the interesting many body aspects in the low temperature regime for 
large values of the local Coulomb repulsion U, whereas leading effects of band-splitting, 
which are present in principle already in the Hubbard-I approximation, and of band- 
deformation can be captured. What is needed for a proper approximation of strong 
correlation effects are summations of infinitely many processes with cumulant vertices 
up to infinite order to cope with the problem of long-time decay of correlations and of 
the infrared problems connected with them. 

The local starting point of the cumulant expansion and its formulation in real space 
open the possibility of a local infinite order resummation, which would be more difficult 
to recognize after unrestricted Fourier-transformation. In the latter fc-space version 
of cumulant perturbation theory, due to the unrestricted site-summations, processes 
taking place on the same site in a diagram, cannot be identified anymore and may 
appear completely uncorrelated. 

Moreover, also individual parts of the partition function in the original denominator 
connected to the same site and necessary for a proper local normalization are hidden in 
the cumulant and cannot readily be identified. A proper conserving approximation of 
infinite order should unite such pieces of local processes in a consistent way. This 
furnishes the guideline for a "locally complete approximation": Collect all those 
diagrammatic contributions to the 1-P-Greensfunction of cumulant perturbation theory 
in real space, where both external legs belong to the same site, not regarding whether 
they belong to the same or to different vertices situated at this site. 

The formulation of a Dyson-equation for the 1-particle-Greensfunction in cumulant 
perturbation theory is straightforward, but needs some care in distinguishing local and 
nonlocal parts of the propagation process: certain intersite-transfers have to be made 
explicit. In a straightforward way one-particle-irreducible pieces can be identified, which 
either are linked to external sources or to transfers. It is useful to amputate local 
factors Q^\z) (the Greens functions of isolated local subsystems) at their two ends, thus 
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Figure 17: Different forms of the Dyson-equation, for the case of one-particle transfer 
only, in cumulant perturbation theory. Care has to be taken in specifying nonlocal 
connections, which is done here via the broken transfer lines made explicit. Single 
unbroken lines relate to the propagators in Free Theory, the hatched rectangular box is 
one-particle irreducible in cumulant (and not in regular) perturbation theory. 



(a) 





Figure 18: Contributions neglected in the selfenergy-matrix in the "locally complete" 
(LC) approximation, one-particle transfer only. Part (a) shows a process, in which two 
loops, based at the same site, are correlated at a different site by a 2-particle cumulant 
vertex. Part (b) shows a process, in which the external particle is extracted at a site 
different from where it was injected. 



defining irreducible cumulant selfenergies ]^ mp \z)- Since the two ends can be situated 
at two different sites, they will generally be fc-dependent after Fourier-transformation. 
Two equivalent forms of the Dyson-equation are visualized in figure 17; their algebraic 
form is: 

Ok(z) = af r) (z) + of T) (z)sir p) W£ (0) (*) + £f T) mp) (^ (0) 

(20) 



= G k (z) + G k (z)t k G k (z) => O k (z)= Chizy'-tk 



with 



1 -1 



(21) 

Ok(z) =a (0) (z)+a {0) (z)^ mp \z)a (0) (z) , (22) 

where E(°)(z) is known from the solution of the local subsystem. The irreducible 
cumulant selfenergy S^ m ^(z) defined above must not be confused with the standard 
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selfenergy 2^ (z) defined via 

Gk(z)=[z-^k\z)-h\~ 1 • (23) 
The connection between them, which can be expressed as 

} (z) = (z) + ^ mp) (z) [l + a (0) (z)^ mp) (z)] , (24) 

sheds some light on the possible momentum dependence of the selfenergy. We expect 
that this may contribute to the question of how to incorporate nonlocal correlations into 
lattice theories (see [27, 60] and refercences therein). 

In a local approximation, the two external links of J^ mp \z) are restricted to be 

situated at the same site, which eliminates the fc-dependence from Y^ amp \z). This type 
of approximation is quite in the tradition of the early effective-site theories mentioned 
above. It brings formal advantages but involves shortcomings concerning the neglect of 
certain nonlocal correlations. 

Still it does not lead to an easily tractable calculational scheme, which is due to 
classes of remaining nonlocal correlations in S^ am ^(z) as indicated in figure 18(a); part 

(b) of this figure on the other hand shows correlations not included in the local form of 
Y^ amp \z). Eliminating diagrams like the one shown in figure 18(a) leads to a restricted 

form Y^ lc \z) of S^ am ^(z), which is at least "locally complete" and may be characterized 
as follows: In a diagram contributing to YS 1c \z) cumulant vertices of order n > 2 (four 
or more external links) are forbidden to appear, since they would correlate two or more 
1-particle loops based at the same site. Therefore 1-particle loops generally will contain 
insertions with two external links (and no more), which are of the same class as those 
contributing to T^ lc \z)] in this way a complete hierarchical structure of independent 
loops remains. 

Amazingly, the locally complete approximation thus defined can quite generally be 
brought into a form suited for a straightforward solution, in principle without further 
approximations or restrictions. It turns out to be equivalent to general formulations of 
the XNCA- or the DMFT-methods [20, 54]. The key to a simpler formulation of the 
locally complete approximation of cumulant perturbation theory lies in the reduction 
of all nonlocal topological elements in a diagram to independent one-particle loops as 
described above. This makes it possible to trace back the fc-summed form of (20), i.e. 

= jf E M*) = jf E - fc] 1 > ( 25 ) 

k k 

to the irreducible part T(z) of the (unrestricted) loop-propagator 

£(*) = ^ZX(*) ' M*)=t ls G k (z)t k (26) 

k 

It turns out that Q{z), as the local one-particle-Greensfunction (i.e. both external 
sources are at the same site), can be constructed as a functional Q[z; Qf(w)] of the 
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irreducible loop-spectrum Qf(oo) = — i/mT(cj+i5), and Q(Q,T) has a simple expression 
via Q and T. Therefore (26) reduces to an implicit equation for T(z) and consequently 

also furnishes solutions for Q(z) , Qk( z ) an d Q(z)- 

The original loop propagator T(z) embodies the possibility that the loop connects 
to its basic site at several intermediate instances due to the unrestricted site-summations 
contributing to G_k(z) in cumulant perturbation theory. The whole loop may be viewed 
as a repetition of irreducible pieces, i.e. loops without such intermediate connections to 
its basic site. The relation between both loop propagators takes the form of a Dyson- 
equation 

T(z) = T{z) + T{z)a{z)T{z) Q{z) = T^)" 1 - T(^)" 1 , (27) 

which was first recognized in [57]. (27) generalizes earlier attempts for implementing 
a "site-exclusion principle" for the propagation of quasiparticles under the influence of 
strong local correlations [3] . 

It is easy to visualize the effect of loops on the dynamics of the basic site of £ amp (z), 
i.e. the one bearing the two external links, when e.g. the picture suggested by direct 
perturbation theory is used. Apart from the fact, that a projection of Qf(uS) onto a local 

1-P-state is used in all formulas instead of the spectrum q^{uS) of band electrons, the 
way to calculate Q^z) is the same as in the impurity problem considered in section 2. 

Insofar, an effective impurity is constructed, and the spectrum of the irreducible 
loop-propagator may be viewed as a matrix of real, frequency-dependent effective 
external fields, constituting a "bath" or a "dynamical mean field". Conclusions about 
Q(z) are to be traced back to the definition Q(z) = G^\z) + AQ(z), AQ(z) = 

Q^\z)]^ mp (z)Q^\z) in (20), reproduced here for the locally complete approximation. 

At first glance it may seem as if the diagrams contributing to Q(z) in cumulant 
perturbation theory via S amp (z) would just reproduce the original contributions of the 
effective impurity problem. However, this cannot be true, since the introduction of 
cumulants also serves the purpose of removing restrictions from site-summations, giving 
Qk(z) the simple form used in (25). QL(z) thus contains compensation terms for loop- 
contributions produced by tk i n the denominator, i.e. in the corresponding geometric 

series with local parts Q(z) and links tk. One possible way of uncovering the relation 

between Q{z) and Q{z) consists in formulating the difference between both quantities 
just as the contribution for loops to be compensated, i.e. $ 

Q(z) - Q(z) = -Q(z)T(z)Q(z). (28) 
J An alternative statement of this relation is 

aw*) = £(*)£(*) . 

which is realized by analyzing the expansions of QL(z) and in terms of irreducible loops T{z). 
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When (27) is inserted into (28) to eliminate Tj^z) in favour of T(z), one obtains 



Q(z)= Qiz^+Tiz) 



a(&T)(z), 



(29) 



thus completing the reduction of (25) to an implicit equation for £ as envisaged above. 

It should finally be remarked that the contributions in cumulant perturbation 
theory to the quantities considered here, which may be classified as connected pieces 
not linked to the external sources in the original picture, are absorbed in the proper 
normalization of spectra and 1-particle-Greensfunctions; they originate from a division 
by the partition function as explained above. Using a consistent locally complete 
summation of these contributions, they become absorbed in the partition function of 
the effective site problem. In direct perturbation theory, for example, this is taken into 
account by properly normalized defect propagators. 

4-2. Periodic Anderson Model: Part I 

The formal development outlined above leads to a result, which constitutes a matrix 
generalization of DMFT. The original concern of LNCA and XNCA was the physics of 
the Anderson lattice model, 



and it should shortly be explained, how a scalar form of the equations is achieved for this 
case. The local subsystem here involves for the simplest case a basis of four one-particle 
states, the two /-states with spin up and down subject to the Coulomb-repulsion U and 
two c-states which do not interact with each other or with the /-states. All matrices 
considered above are four by four, with £ being spin-diagonal and transferring only c- 
and /-electrons to nearest-neighbour c-states with respective matrix elements t and V. 

Although unphysical in most cases, a purely local hybridization V is often used 
for simplicity; it can be treated in close analogy to the nearest-neighbour case for the 
reason explained in the following. Since the c-electrons remain noninteracting, Wicks 
theorem can be used for them. Consequently, no cumulant vertices of order n > 2 exist 
with links to c-Greensfunctions. Q^\z) is block- diagonal with a diagonal c-block and 

a /-block, and likewise is which is used in the amputation of vertices. 

The block-structure mentioned is 2 x 2 with respect to spin; without magnetic 
splitting and with spin-preserving transfer and hybridization one may fix a spin-direction 
a and treat these blocks as scalars. As a consequence, S amp (z) and also AQ^z) have 

nonzero matrix elements only in the diagonal, i.e. only Gff a and G cccr (z) = G® C(J (z) 
enter the calculation according to the definition in (20). Nondiagonal elements come 
into play only via T(z), since propagation along a loop can mix c- and /-states. (29) 
now gives after matrix-inversions: 




Tffa(z) = G ffa (z) 



-1 




cc<j\Z) ^ ~cc(j\Z 



-1 
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T cfa (z) = T fca (z) = , N a {z) = G ffa (z)G cca (z) ~ G efa (z)\31) 

These matrix elements all would have to be used if the four states were locally correlated. 
However, in the simple form of the Anderson model the local c-states do not directly 
influence the /-state dynamics. Therefore one can reduce the local problem to one 
of /-states by combining site-irreducible loops in a way that only a restricted form of 
irreducibility with respect to the /-states on the basic site is realized, i.e. the c-level on 
this site is treated like a different site. 

The combined loop-propagator connects local /-states and obeys a generalized 
Dyson-equation of the form fjjf\z) = f ff + f fc G cc f cf + f fc G cc f cc G cc f cf + ... and 
hence may be summed up to: 

ffa ^ ' /M ' l-G cca {z)T cca {z) K ) 

If now the quantities T on the r.h.s. are replaced via (31), one obtains the scalar 
equivalent to (29): 

Tf f 1\z) = Gff^z)- 1 - GffAz)- 1 ■ (33) 

As is clear from this construction, the spectrum of T^ re ^ is to be used in the effective 
site problem, i.e. Gff a (z) = Gff a z\ Qf(red)(uj) . 

Also the inverse of the last term in (33), i.e. Gff a (z), has to be calculated 
via (25) using a matrix-inversion, which involves the non diagonal transfer-matrix 
tk = tkX + Vkg^. After explicitly formulating this step the reduction from a matrix- to 
a scalar form of the theory for the Anderson-lattice model is completed. 

Whereas there might exist easier ways of setting up the XNCA/DMFT- 
self-consistency cycle for this model [20, 62], the above argumentation generally 
demonstrates the connection between the universal matrix-formulation and possible 
reduced schemes. If, for example, electrons in c-states would interact locally (but again 
not with those on /-states) , the matrix problem would reduce to two scalar problems of 
type (28), which would be coupled only via the lattice-summation, i.e. via the fc-sums in 
(25)). This also points to a possible treatment of more general models with inequivalent 
sites. 

Finally, since formally the local c-state acts like a different site on its /-states, a 
local hybridization can be treated in the same way as outlined above: The first and 
last transfer step, represented in (26) by factors now carries a fc-independent matrix 
element V. This only enters the calculation in a modified = tk 1 + Vg^ to be used in 
(25). 

The Anderson-lattice model furnishes a good testing ground for the quality of 
impurity-solvers because of the particular impact of coherence in the half-filled case. 
In the symmetric situation of the simple version without orbital degeneracy considered 
above, with two electrons per site, the Luttinger-theorem predicts a Fermi surface filling 
the whole first Brillouin zone, which should lead to the formation of an excitation 
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gap due to Bragg scattering in the quasiparticle-DOS for temperature T approaching 
zero [59, 61]. This signature of onsetting coherence is hard to recover in approximations, 
since it requires a pronounced structure in the selfenergy. 

On the one hand the increasing lifetime of quasiparticles with T — > near the 
Fermi level requires the (near) cancellation of the term (see last section) in the 
local selfenergy by the builtup of scattering during propagation along loops through the 
lattice (see above), and on the other hand formation of a gap should go along with a 
narrow peak in ImS signalling strongly increased resonant scattering. This implies a 
fine balance between different contributions, which is easily destroyed by inconsistent 
approximations. 

4.3. Hubbard Model 

In order to elucidate the effect of increasing lifetime as T —> one may consider the 
Hubbard model, 



where at half-filling and zero temperature scattering should be absent near uo = in 
the Fermi-liquid phase, but which should not develop a coherence gap: Since according 
to Luttingers theorem the Fermi surface now lies well inside the first Brillouin zone no 
Bragg scattering should be effective there. 

Figure 19(a) shows the corresponding local DOS, obtained within different 
approximations, two of them within the locally complete scheme (DMFT) using SNCA 
and ENCA as impurity solvers, respectively. 

A reduction process from the matrix-formalism analog to the Anderson lattice 
model is not necessary for the Hubbard model, since the latter is of scalar type from 
the outset. 

The other approximations are Hartree-Fock, effectively meaning U = in the half- 
filled case, and Hubbard-I ("Free Theory"). The last two cases show no temperature 
dependence, whereas (3 = 100 and f3 — 10 are chosen for SNCA and ENCA in order 
to produce comparable heights of the many-body resonance, which should reach the 
Hartree-Fock value at T = 0. 

In a fc-resolution this DOS produces the quasiparticle-bandstructure, as shown 
for the ENCA-calculation in figure 19(b) along the [11 Indirection of the simple cubic 
Brillouin zone. At high excitation energies the two split bands, which in Hubbard- 
I approximation contain sharp resonances with reduced spectral weight, are so much 
washed out that a concept of band electrons can hardly be justified here; this is in 
accord with the conclusion in the last section about local scattering near an impurity. 
The narrow band of pronounced QP-resonances around uj = ji = shows no splitting 
and gap- formation. The corresponding decrease of scattering is shown in figure 19(c): 
In lowering the temperature from /3 = 2to/3=10 the imaginary part of the self-energy 
forms a steep and nearly quadratic minimum as a sign of Fermi-liquid formation. 




(34) 
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Figure 19: One- particle excitation spectrum of the Hubbard model, calculated within 
XNCA/DMFT, using ENCA as impurity solver, for the half-filled case with nearest- 
neigh bour hopping in a 3d-sc lattice. Parameter values e = —2, £7 = 4 and (3 not too 
large favour a metallic phase with a Fermi liquid. In part (a) the local DOS is compared 
for different approximation schemes, also including Hartree Fock (resulting in a tight- 
binding band of width 6) and the Hubbard-I approximation. Beside ENCA with j3 = 10 
also SNCA with (3 = 100 is used as impurity solver. Part (b) shows the fc-resolved 
excitation spectrum along the [11 Indirection, from which a quasiparticle bandstructure 
may be derived, within XNCA/DMFT-ENCA. The imaginary part of the selfenergy 
(absolute value) in part (c) visualizes the formation of the Fermi liquid with decreasing 
temperature (fcs)T = f3~ 1 . 
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Figure 20: One-particle excitation spectra of band (o) electrons and of local (/- 
) electrons for an Anderson-lattice model, calculated within XNCA, using CA1 as 
impurity solver, for the half-filled case (two electrons per site), with a tight-binding 
c-band of width 6 in a 3d-sc lattice and local hybridization. Parameters are = 
-1.5, U = 3, f3 = 15 and A A = ttV 2 q { $(0) = 0.3. Parts (a) and (b) show the local 
density of c-and f-electrons, respectively, both in comparison with a Hartree-Fock result 
and a calculation within "Free Theory". Part (c) contains the imaginary part of the 
"local selfenergies" (absolute value) for c-and f-electrons. The spikes seen in the middle 
of the gap region can cause numerical problems. 
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Figure 21: The fc-resolved excitation spectra of c-and f-electrons, respectively, for the 
same parameters as in 20 are drawn along the [11 Indirection. They demonstrate, like 
parts 20(a) and (b), the formation of hybridization pseudogaps, smeared by lifetime 
effects, in the high-energy region and of a narrow and complete coherence gap at the 
Fermi energy. 
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Quite generally it can be stated, that the ENCA impurity solver for the Hubbard 
model produces good results and accomplishes an acceptable tradeoff between accuracy 
and numerical effort. For not too low temperatures it thus represents a reliable and 
usable alternative to impurity solvers like NRG, QMC, MPT, ...known from the 
literature (see for example [63, 22]). 

4-4- Periodic Anderson Model: Part II 

The Anderson-lattice model behaves differently than the Hubbard model, as shown in 
figure 20 and 21. Parts 20(a) and (b) contain the local one-particle excitation spectra for 
c- and /-electrons, respectively, each calculated with three approximations of increasing 
complexity. The tight-binding approach for the band (c-) states treats all interactions 
on the mean-field level and furnishes the connected curve known from the Hubbard 
model with its edge-like van-Hove-singularities at uj = ±1 and square-root band edges 
at uj = ±3. 

Hybridization with the correlated local /-states at one-particle energies uj = = — 1 
and uj = €£ + U = +1 produces the gaps visible in the result of the Free Theory, 
which furnishes three disconnected bands with variable spectral weights. The gaps 
are somewhat displaced by level repulsion, which is an effect of hybridization, too. 
Interactions, much better taken into account in the locally complete approximation 
using a CAl-impurity solver, wash out the two gaps and produce a repulsion of c-weight 
away from the Fermi level uj = as a consequence of the formation of the many-body 
resonance with predominant /-character. This resonance is clearly seen in figure 20(b), 
where also the two local one-particle levels of the isolated /-shell and the spectrum of 
the Free Theory are shown. In the latter, two gaps are recognized as counterparts of 
those in figure 20(a); the /-states acquire dispersion through mixture with the band and 
share, in corresponding regions, the effect of level-repulsion and gap-formation. 

Interestingly, the CAl-impurity solver at the low temperature considered, i.e. 
(3 = 15, is able to describe the formation of the gap in the narrow region of quasiparticle 
states near uj = 0, which was to be expected as a consequence of coherence in the 
Anderson-lattice. This effect is connected with a strong increase of scattering at uj = 0, 
see figure 20(c) where imaginary parts of the "local" self energies 

^aa(z) = z - e a - (Jy G aa ka(z^j , (z = Huj + iS , a = c, /) (35) 

k 

are shown; the corresponding very narrow spikes at uj = can cause numerical problems. 
In approximate impurity solvers usually convergence problems of the XNCA/DMFT- 
cycle are observed. Whereas the original gaps become smeared also in the /-spectrum, 
the coherence gap should become perfect in the limit T —> 0. This leads to the 
narrow gapped quasiparticle band structure, to be seen in the fc-resolved spectra of 
figure 21(a) and (b), and to the broad structures, smeared out by the interactions, at 
higher excitation energies. 
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Finally it has to be mentioned that stability of the Fermi-liquid phase might not 
be thermodynamically stable for all parameter values used in our calculations for the 
Hubbard- and the Anderson-lattice model. Susceptibilities can point to instabilities 
towards other possible ground states [65], which is among a variety of methods being 
applied to the investigation of phase diagrams. 

5. Conclusion and outlook 

The foregoing sections have demonstrated the considerable progress, which has been 
made in the development of impurity solvers via direct perturbation theory and their 
application to impurity- and lattice-problems with strongly correlated electrons. Our 
presentation has emphasized an unified view on several approximations of this kind, 
which have been proposed in the past, and on a new one, the CA1, discussed here for 
the first time. 

All of these approximations can be characterized in a systematic fashion as 
skeleton expansions in terms of time-ordered local perturbational processes along the 
lines laid out by [31]. As such, they furnish coupled implicit integral equations for 
propagators, which in general have to be solved numerically; this gives rise to the 
notion " semianalytic" . ENCA, SUNCA, FNCA and CA1 include different classes of 
vertex corrections; the first three of these approximations reduce to SNCA, the version 
of the old NCA without any vertex correction applied to finite U. CA1, on the other 
hand, contains fully crossing vertex corrections of fourth order in the hybridization. 

It turns out that for the quality of the approximation it is important to include 
ladders for repeated particle scattering and higher oder vertex corrections in a well- 
balanced way. This is apparently accomplished best by the CA1, which however does 
not iterate special subclasses to infinite order like SUNCA and, more generally, FNCA. 
Comparison with NRG-calculations in the spirit of Wilsons approach, reveals that even 
CA1 has deficiencies at low temperatures and excitation energies. At higher energies, 
however, the situations is reversed: In an overall view taking into account the complete 
spectral region, the semianalytical impurity solvers, with the possible exception of 
SNCA, perform quite well, and even the ENCA, as the least complicated of them, 
may be used for qualitative investigations. 

The unified view developed here also concerns the construction of approximations 
for lattice problems and the use of the impurity solvers therein. It was shown that 
appropriate choices of local building blocks, each containing a set of internally correlated 
one-particle states, and a selection of paths for propagation through the lattice can be 
consistently combined in a matrix-formulation for a calculation of partition function 
and Greensfunctions. 

It was explained how in certain simple situations, as e.g. encountered in the 
Hubbard model or the Anderson-lattice model with noninteracting band states, the 
formalism reduces to a scalar one and how a locally complete selections of local processes 
then leads to the well known XNCA-and DMFT-approximations. 
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Our starting point was a cumulant expansion for all local n-particle vertices, which 
renders the application of the linked cluster theorem and of unrestricted site summations 
possible, thus enabling a convenient fc-space representation of quantities. Generalized 
dynamical fields have been introduced and traced back to matrix-propagators along 
closed loops. The neglect of all cumulant vertices of order n > 2 leaves as a natural 
"Free Theory" the matrix-generalization of the Hubbard-I approximation, which can 
readily be evaluated in explicit form. Restricted selection of cumulant vertices to 
finite order allows to define e.g. Hartree-type approximations [54, 56], which are not 
yet well investigated but supposedly are of restricted usefulness for the regime of low 
temperatures and excitation energies. 

The locally complete approximations, on the other hand, combine advantages 
of real space as well as of fc-space formulations with a better treatment of infrared 
divergencies through infinite order. Nowadays this can be implemented by a variety of 
local impurity solvers, among which we have concentrated here on the class based on 
direct perturbation theory. As applications of the formalism we have presented local 
and fc-resolved one-particle-excitation spectra for Hubbard- and Anderson-lattice model 
and have discussed characteristic similarities and differences. It has proven useful to 
connect this discussion with the foregoing treatment of the SIAM as the prototypical 
effective impurity. In particular, the formation of a Fermi liquid could be illuminated 
in this way, emphazising a local point of view. 

We will conclude with a short perspective on possible future developments on the 
basis of our local approach and with some critical remarks about its shortcomings. 
Improvements of semianalytical impurity solvers could be based on CA1 and proceed 
along directions laid down in the simpler case of the U = oc-version of SIAM [41, 44]. 

Two different approaches could be combined to develop such a CA2-theory. The 
fully crossing fourth order vertex corrections could be reinforced by certain diagrams 
of sixth order like in figure 4(b), and possibly iterated further, which had proven as 
beneficial for spectral properties at large U in the frame of the PNCA [41]. The 
CTMA [44, 46], on the other hand as the second of these theories for U = oc, stresses 
the role of long ladders of crossing particle lines. 

A generalization to finite U can be accomplished by solving the system of four 
coupled T-matrix equations shown in figure 22. The resulting T-matrices then are to 
be inserted into the fully crossing vertex corrections contained in CA1; they replace 
parts shown as lowest order contributions envisaged in figure 3. The first T-matrix of 
figure 22 e.g. replaces the diagram in figure 4(a), leading to the sequence of figure 5(b). 
Both measures together should again be well balanced in the sense discussed above. 
Although we expect further essential quantitative improvements by such additions to 
CA1, the resulting CA2-impurity solver can still not be expected to be perfect. 

It lies at the heart of the infrared problems in SIAM or its generalizations that no 
treatment based on a restricted selection of perturbational processes to infinite order can 
adequately describe the complete energy range down to uj = 0. From a practical point 
of view, however, for many purposes this will not really be necessary, for example when 




Figure 22: Higher order corrections to CA1 can be implemented via the four T-matrix 
equations shown graphically. The first of these was considered in figure 5. The other 
three equations would generalize U = oc-theories like PNCA and CTMA to finite values 
of U (CA2-project). 

other types of correlations in concentrated systems intervene, e.g. producing magnetic 
states. One problem will be left anyway, even in this case: The numerical expense 
connected with approximations like CA1, and even more so with a hypothetical CA2, 
is considerable. Certainly there will be a need for improved algorithms or for the 
use of parallel computing. At the end, these higher semianalytical impurity solvers 
may turn out impractical compared with e.g. Quantum-Monte-Carlo or NRG-methods. 
Their ability, on the other hand, to describe the regime at higher temperatures and 
excitation energies very well, is already shared by the ENCA, which does not need so 
much numerical effort. 

Impurity solvers in general serve as a key ingredient in the local approach to 
lattice systems. Current research aims at the inclusion of more realistic local building 
blocks containing several orbitals or clusters of sites and at a better treatment of 
nonlocal correlations. With an optimal selection of a localized basis of Wannier-states 
screened direct interactions are hopefully short ranged and can either be treated within 
a small cluster as a local building block or via an extension of the perturbation in the 
Hamiltonian to double-transfer between neighbours, as indicated in section 2. 
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This, as well as a multi-orbital situation, can in principle be handled by the matrix- 
propagator formalism presented in this paper; it has been demonstrated in section 4 
how matrix-loop propagators systematically serve to define appropriate dynamical fields 
acting on effective sites. Long ranged correlations, however, need an extension of the 
locally complete approximation for their proper incorporation. 

The cumulant approach outlined above, can serve as a basis for such generalizations. 
The nonlocal processes additionally to be taken into account in e.g. selfenergies or 
susceptibilities connect local cumulants with four external legs or more on different sites 
via more than one loop. The resulting correlations between propagating particles or 
holes and between the states of local blocks on different sites may generate different 
types of long-range order, e.g. phases of itinerant or local-moment magnetism. 

Technically, the difficulty in setting up such appropriate nonlocal theories seems to 
lie in a consistent choice of diagram classes as skeletons without overcounting, so that 
basic requirements are fulfilled, such as conservation laws and analytic properties like 
positivity of excitation spectra. This requires, apart from a strict use of cumulant 
subtractions, the consequential implementation of the concept of irreducibility [21]. 
Hopefully, work along these lines will soon lead to improved forms of e.g. band structure 
theories for correlated electrons and will shed more light on the mechanisms behind the 
formation of exotic ground states in transition metal compounds. 
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